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The asymptotics of an eigenfunction-correlation 
determinant for Dirac-(5 perturbations 

Martin Gebert 


Abstract. We prove the exact asymptotics of the scalar product of the ground 
states of two non-interacting Fermi gases confined to a 3-dimensional ball Bl 
of radius L in the thermodynamic limit, where the underlying one-particle op¬ 
erators differ by a Dirac-5 perturbation. More precisely, we show the algebraic 
decay of the correlation determinant \ det ^ ^ 

as M,L —>■ oo and N/\Bl\ p > 0, where ipf and denote the lowest- 
energy eigenfunctions of the finite-volume one-particle Schrodinger operators. 
The decay exponent is given in terms of the s-wave scattering phase shift 
({E) := 5^{y/E)l7v^. For an attractive Dirac-5 perturbation we conclude that the 
decay exponent ^|| arcsin |r(i?)/2| ||hs found in |GKM014!| does not provide a 
sharp upper bound on the decay of the correlation determinant. 


1. Introduction 

We consider the asymptotics of the scalar product of the ground states of two 
non-interacting finite-volume iV-particle Schrodinger operators in the thermody¬ 
namic limit approaching the particle density p{E) > 0 corresponding to the Fermi 
energy E > 0. Here, the underlying one-particle Schrodinger operators are the neg¬ 
ative Laplacian in 3-dimensional Euclidean space and the negative Laplacian with a 
Dirac-(5 or zero-range perturbation sitting at the origin. We restrict this pair to the 
ball Bl(0) of radius L and are interested in the L-asymptotics of the scalar product 
of the ground states of the corresponding two non-interacting Wparticle operators, 
which we call ground-state overlap in the sequel. Using the representation of the 
ground states as Slater determinants, we see that the ground-state overlap is the 
following correlation determinant 

( 1 . 1 ) 

In this note, we are interested in its thermodynamic limit, i.e. increasing L and 
N G N simultaneously such that N/\Bl(0)\ —)• p{E) > 0, where p{E) denotes 
the integrated density of states of the negative Laplacian at the energy E > 0. 
Here, ipj and are the eigenfunctions belonging to the N lowest eigenvalues of 
the restricted operators, which we call Ef^ and and (•, •) denotes the scalar 
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product in L^(i?i(0)). Anderson claimed in [And67] that in the case of a Dirac-5- 
perturbation the determinant admits the asymptotics 

( 1 . 2 ) 

as iV,L —>■ oo, N/\Bl{0)\ —>■ p{E) > 0, where C{E) ■= and 6 refers to the 

s-wave scattering phase shift. This algebraic decay of the ground-state overlap is 
called Anderson’s orthogonality catastrophe in the physics literature and we refer 
to |GKM14] for further references. 

The starting point of the proofs of previous rigorous results is the following 
expansion of the determinant 

OO . 

hf f = - E. (13) 

n=l 


valid for appropriate choices of N, where and denote the A^th and 

(A -|- l)th eigenvalue of the finite-volume operators Hl and Ha,L, see [GKM014J . 
Thus, estimates on the correlation determinant are closely related to asymp¬ 
totics of products of spectral projections given in (11.31) . Gonsidering only the n = 1 
term in m, the hrst rigorous bounds on were proved in [KQS13] valid 
for 1-dimensional systems and short-range perturbations. They found the upper 
bound on < L~^ with the decay exponent ^{E) := ^||T(A)/2||^g, where 

T refers to the scattering T-matrix of the corresponding infinite-volume opera¬ 
tors, and a non-optimal lower bound. Later in |GKM14] the same upper bound 
j{E) := ;j^||r(A)/2||^g was deduced for quite general pairs of Schrodinger oper¬ 
ators in arbitrary dimension, which differ by a sign-definite potential. Taking all 
summands in (jl.3p into account, [GKM014] proved an upper bound with the de¬ 
cay exponent 'y{E) := ^||arcsin |r(A)/2| ||gg in the general setting discussed in 
[GKM14] . Let us point out that these previous results concern upper bounds and 
are also valid for special choices of thermodynamic limits only. Here, in the toy 
model of a Dirac-h perturbation we provide the exact asymptotics of the correla¬ 
tion determinant and we consider arbitrary thermodynamic limits approaching a 
particle density p > 0, see Theorem 12.II below. We show this using a representation 
of the ground-state overlap other than m, which is valid for rank-l-perturbations, 
i.e. 




A A IfEdTl-d 


(1.4) 


where and pj are the eigenvalues of the pair of the finite-volume Schrodinger 
operators, see Section [3l This formula is known in physics literature and goes back 
at least to [T085] . Using the latter formula, we give a straightforward proof of the 
algebraic decay with the exponent C{E) = as Anderson predicted. 

It turns out that the decay exponent is equal to the one found in [GKMOI4] in the 
case of a repulsive Dirac-h perturbation only, i.e. C{E) = ^||arcsin |T(A)/2|||^g. On 
the other hand, we obtain C{E) > ;j^||arcsin |r(A)/2|||gg, for an attractive Dirac-(5 
see Remark 12.31 below. Hence, the decay exponent ^||arcsin |r(A)/2| ||gg does not 
provide the exact asymptotics of m- 

Recently, [KOS15] proved the asymptotics of a shifted correlation determinant 
for one-dimensional models with a perturbation by a magnetic field. A related 
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problem, which we mention for completeness, is considering the asymptotics of 
products of spectral projections of infinite-volume operators, similar to (II.3p . This 
was done in the proof of |GKM014| and extended in [FP_14j . 

2. Model and results 

We start with the operator — Aq : —)• which has deficiency 

indices (1,1). Therefore, —Aq gives rise to a one-parameter family of self-adjoint 
extensions which we index by a G R and denote by — Aq, see [AGHHOSi Chapter 1]. 
We refer to — Aq as the negative Laplacian with a Dirac-(5 perturbation sitting at the 
origin 0 of strength a. Throughout, we consider for a G R the pair of Schrodinger 
operators 

H := —A and '■= —Aq (2T) 

on the Hilbert space Ti = L^(R^), where —A is the negative Laplacian. More 
precisely, following [AGHHOSi Chapter 1], the operators H and Ha admit a decom¬ 
position with respect to angular momentum. Thus, there exists a unitary U such 
that both operators transform into the direct sum 

UHU* = 0 0 and UHaU* =0 0 (2.2) 

-^ENo -^ENo 

where : L^((0,oo)) D dom(/i^^^) —>• L^((0,oo)) and for ^ ^ 1. In 

the i = 0 case the operators are given by 

dom(/i°) = {/ G L^((0,oo)) : /,/' G AC'ioc((0, oo)); (2.3) 

/(0+)=0; f”GL\{0,oo))} 

dom(/iQ) = {/G L^((0,oo)) : G ACioc((0,oo)); (2.4) 

— 47ra/(0-|-) -|- /^(0-|-) = 0; f” G T^((0, oo))}, 

where we denote by AC'ioc((0, oo)) the set of all locally absolutely continuous func¬ 
tions. Thus, the difference of H and Hq. takes place in the lowest angular momentum 
channel via a different boundary condition at 0 which we parametrise by a G R. 
In the following we are interested in the restrictions of these operators to the ball 
Bl{0) of radius L around the origin 

Hl ■■= -Al and Ha,L ■= -^a,L- (2.5) 

Here, — A^, is the negative Dirichlet Laplacian on Bl(0). The operator —A^^l cor¬ 
responds to the restriction of the operator — Aq imposing Dirichlet boundary con¬ 
dition at L in each angular momentum channel, i.e. also the restriction of —Aq 
to Bl{0) with Dirichlet boundary conditions. Thus, Hl and Ha,L differ as well 
as before in the lowest angular momentum channel only by a different boundary 
condition at 0. We call the corresponding operators in the £ = 0 channel, i.e the 
restrictions of and to the interval (0, L) with Dirichlet boundary condition at 

L, 

hi and h^L- (2-6) 
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Using standard results for regular Sturm-Liouville operators, we obtain for all z G 
p{h^) n p{h^ j^) a vector ^ G L‘^{Bl{0)) such that the resolvents satisfy 


1 1 






a 

' L , z \ 


(2.7) 


Thus, is a rank-1-perturbation of in the resolvent, and the same is true 
for the pair and Hi. We point out that the perturbation is not compactly 

supported since ^ is L dependent. Moreover, the compactness of the resolvents 
oi Hl and H^ l imply that both Hl and Ha^L have discrete spectra. We write 


Af ^ Af ^ • and /if ^ /i^ ^ • 


( 2 . 8 ) 


for their non-decreasing sequences of eigenvalues, counting multiplicities, and 
((/jf )jgis} and {'4>k)ke.f^ for the corresponding sequences of normalised eigenfunctions, 
where we choose the same eigenvectors for Hi and in any angular momentum 
channel i ^ 1. This choice ensures that the eigenfunctions of Hi and differ 
in the lowest angular momentum channel only. Let us point out that in the case of 
a < 0 there exists precisely one negative eigenvalue pi = — (dvra)^ for the infinite- 
volume operator Ha, respectively see [AGHHOSl Chapter 1]. Dirichlet-Neumann 
bracketing implies ^ /if > where /if ^ denotes the negative Laplacian on 

(L, oo) with Dirichlet boundary condition at L. Thus, in the case of a < 0 we obtain 
the uniform lower bound on the finite-volume operators 

Ha,L ^ — (dvra)^ and equivalently /if^^ ^ — (dvra)^. (2-9) 

Let A G N. In the following we are interested in the correlation determinant 

SH ■■= P-10) 

The main result concerning is the following. 


Theorem 2.1. Let a G 

subject to 


E > 0 and N(^.'^{E) : M_|_ —)• N an arbitrary function 




( 2 . 11 ) 


|Bl( 0)| ' ^ ■ 87r3 ’ 

i.e. p denotes the integrated density of states of the operator —A. Then, the corre¬ 
lation determinant corresponding to the pair Hl and Ha,L admits the asymptotics 


|^IV^(E)|2 ^ ^-^sl(VE)+o{l) 


as L ^ oo, equivalently, 


In L 

and 6a is given by Definition \2.2\ below. 


lim 
L —^ oo 




( 2 . 12 ) 

(2.13) 


Definition 2.2 (Scattering phase shift). Let /c > 0. Then, the scattering phase 
shift is defined by 


6a{k) := 


arctan 

< 



k 


for a ^ 0, 


TT — arctan 


for a ^ 0 


(2.1d) 
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where we use the convention arctan 


TT 


:= — for /c > 0. 


Remarks 2.3. (i) The separate definitions of the phase shift are reminiscent 

to the existence of a negative eigenvalue whenever a < 0 and Levinson’s theorem. 

(ii) Due to the nature of a Dirac-<5 perturbation in 3 dimensions the same result 
is apparently valid for the corresponding problem on the half-axis. 

(iii) We emphasise that we allow arbitrary thermodynamic limits approaching 
the particle density p > 0. 

(iv) The o(l)-error in (j2.12D depends on the particular choice of the thermo¬ 
dynamic limit. To see this, we refer to equations (I4.38p and (|4.39p in the proof of 
Theorem 12.11 In particular, we think that the error cannot be improved allowing 
arbitrary thermodynamic limits. 

(v) In [GKM0141 Theorem 2.2] an upper bound on the ground-state overlap 
is proved for quite general pairs of Schrodinger operators which is valid for subse¬ 
quences only. More precisely, they prove for a subsequence 

,JVi,(E)| 


lim sup ■ 

L —^ cxD 


In Sr 


liiL 


^ - 


l{E) 


where 


7 (L 1 ) := ^||arcsin|r(L;)/2||| 


vr^ 


HS 


(2.15) 


(2.16) 


and T denotes the scattering T-matrix. Since we consider here s-wave scattering, 
we restrict ourselves to the lowest angular momentum channel. In this case, T{E) 
is a complex number and \T{E)/2\ = sin{6a{y/E))- Now, computing 7(L1) yields 


liE) = 




for \6a{VE)\ ^ f 


TT^ 


arcsin ( sin((5a(\/^))) for ISaiVE)] ^ 


(2.17) 


Thus, in general the decay exponent 'y{E) does not provide a sharp upper bound 
on the correlation determinant whenever the phase shift is bigger than 7r/2. In our 
model this is equivalent to a < 0 which we refer to as the attractive case. 


The proof of Theorem 12.11 follows from a different approach than the one made 
in [GKM14] and [GKM014], i.e. we do not use the representation (m in this 
article. Here, the key is the following remarkable product representation of the 
determinant in terms of the eigenvalues of the finite-volume Schrodinger operators. 
To our knowledge, this was hrst stated in [T085j . 


Lemma 2.4. Let N Then, 








(2.18) 


We start with proving this product representation for general pairs of compact 
operators which differ by a rank-1-perturbation in Section [3l We apply this to our 
setting in Section [4] and prove Theorem 12.11 
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3. Representation of the ground-state overlap 


In this section we prove a quite general representation for determinants of eigen¬ 
vectors of pairs of operators which differ by a rank-1-perturbation. The main result 
in this section, Theorem 13.11 will be the key to the proof of Theorem 12.11 

Let Ti. he & separable inhnite-dimensional Hilbert space and A : TL ^ TL he 
a compact, linear and self-adjoint operator. Moreover, we assume A ^ h with 
ker(H) = {0}. We define 

B-.= A + \4>){(t>\ (3.1) 

for some 0 ^ (j) ^ T-L. We write ai ^ 02 ^ • and /?i ^ /32 ^ • for the non¬ 

increasing sequences of eigenvalues of A, respectively B and denote by 
and {'ipk)k£n corresponding normalised eigenvectors. Since A and B differ by a 
rank-l-perturbation, the min-max theorem implies that the eigenvalues interlace. 
We assume the following condition on the eigenvalues 

CO 

^ |/3n - Gnl < 00 . (3.2) 

n=l 

Moreover, for simplicity we also assume the following strict interlacing condition 

Pi > ai > /32 > a2 > ■ ■ ■ . (3.3) 

In particular, Pk 7^ Oij for all j. A; G N. Furthermore, the above implies cyclicity of p. 
Assumption (13.3p is not necessary but simplifies notation and computations. In the 
general case one has to consider the restriction to the cyclic subspace generated by 
the perturbation p. But the application in mind will satisfy the interlacing condition 
(j3.3p . therefore, we assume it. 


Theorem 3.1. Let N gN. We assume eonditions (1^ and (123]) to hold. Then, 


N 


det 


Tj: 




n n 

j=l k=N+l 


\Pk - aj \ \ak - Pj 
\0ik - Oij \ \Pk - Pj 


(3.4) 


Proof of Theorem 13.11 We use the eigenvalue equations and assumption (13.31) 
to obtain for all j,k GN 


{Tj^i’k) 


{ipj,p){p,pk) 

Pk - oij 


Hence, the multi-linearity of the determinant implies 

^ ^j {Tj,P){P,'>Pk) \ 2 

V Pk- aj ) 

N N 

j=i k=i 



det 


1 


Pk - aj 


l^j,k^N 


(3.5) 


2 


(3.6) 














THE ASYMPTOTICS OF AN EIGENFUNCTION-CORRELATION DETERMINANT FOR DIRAC-d PERTURBATIONS 


Now, the remaining determinant can be computed explicitly. We use the Cauchy 
determinant formula to evaluate this, see e.g. |Weyl3| Lem. 7.6.A], and end up with 

»=( n n Kv’i. '' 

^i=ifc=i 

Corollary 13.31 below yields 


nj,A;=l \/^k Oijl 


(3.7) 


N oo 


62»= nn 


\ai - Pk\ 




N oo 

nn 


— aj 


l^k 
N oo 

n n 




N 

n 

j,k=l 


\Pk - Pj\ - «fc| 
\Pk - ajP 


|/3fc - aj\ \ak - /3j\ 


(3.8) 


\Pk — PiWdi — CKkl 
j=lk=N+l r'Jll J 

This gives the claim, where we remark that by assumption (13.2p all products in the 
latter converge absolutely. □ 

To complete the proof, we continue with computing the resolvents of the oper¬ 
ators A and B in terms of their eigenvalues. 

Lemma 3.2. We assume (j3.2p and (13.31) . Then, there exist a, 6 G M with ab = —1 
such that 


(i) for all z G q{A) 


/A ^ A\ I ^ YJ Pk — Z 


k=l 


(ii) for all z G q{B) 


{ 4 >, 


1 


y>-i = i.n 


ak- z 


an - z 


B-z^' ■ 

Corollary 3.3. Let j,k G N. Under the assumption (13. 2p and (13.31) 


\{Tj,P){i^k, 4 >)\'^ = \Pj - ajl \ak - / 5 fc| ( n 


— a, 


\ai - Pk\ 
\ai -aj\J y \/3i - I3k\ 

l^j Ij^yk 


n 


(3.9) 


(3.10) 


(3.11) 


Proof of Corollary 13.31 Using Lemma 13.21 we compute the residue of the re¬ 
solvents 


2 1 
\{Tj^4>)\ = lim {o^j - - 4>) 

z^aj ■' A — Z 


= lim {oj — z) a IT ^^ = a (/?, — 


(«/ - z) 


{Pj-aj)Yl 


iPi - «i) 

fJl im-oij) 


(3.12) 


and along the same line 


\{fjk,(t>)? = b{ak- Pk) n • 


l^k 


(3.13) 
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Taking the absolute value and using |a6| = 1, we get the result. 


□ 


Proof of Lemma 13.21 First note that by assumption (|3.2I) the sequences 


N 

n 


Pk- z 


and 


N 

n 


On - 


(3.14) 


■ k=i ^ ^ V / Nen 

converge locally uniformly for all z E Q{A)r]Q{B), see jKno96i Thm. 252]. Therefore, 
the limits 


Fiz) := n 


On - Z 


and G{z) := 

fc=i 


Pk- z 


(3.15) 


z ■ ' ak- z 

a=l k=l 

are well-defined analytic functions on q{A) H q{B), which fulfil FG = 1. Due to the 
locally uniform convergence, the derivative of F satisfies 

N N 




an — z d ai — z 


N^oo Bn — Z dz Bj — Z 

1=1 n=l 
n^l 

N N r, 

y TT ~ z ai — PI 

^ A^Too ^ M (A - ^)2 

n^l 


N 

F{z) lim (-— 

N^OO V A “ 


1 


IBi- z ai - z 


for all z E ^?(^) n q{B). We apply Lemma 13.41 below and obtain 

Now, the resolvent identity implies for all z E ^?(^) H q{B) 

1 1 




B - z A- z A - z ^B - z 
which provides the equality 

1 , 1 1 , 


•) 


1 - 

Inserting this into (13.171) . we see that F solves the differential equation 


F'{E) = F{E) 


{4’j 


-(A 


B-E 


(t>) 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 


at least for all E E q{A) n q{B) n M. On the other hand the resolvent of B 
is analytic in q{B) and the function t i— {4>, t < 0, solves the 

above ODE (I3.20p as well. Now, the general solution to this ODE is f{t) = 

xoexp^ /tosome initial condition {to,Xo). Note 

that the functions 1 1 —)• F(t) and t i-A (0, —1 are non-zero, thus ((/>, — 1 = 

cF{t) for some c 0. This and the identity theorem for analytic functions give the 
claim. Equation (]3.9D follows from E{z)G{z) = 1 and the identity 

for all z E ^?(^) n q{B) which is a consequence of (13.181) . □ 
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Lemma 3.4. Let z E H q{B). Assume (13.3p . Then, we obtain the following 
identity 



(3.22) 


Let us point out that in the finite-dimensional case the above equality follows 
directly from the resolvent equation, (|3.18p . Nevertheless, the infinite-dimensional 
case is a bit more involved due to convergence issues. 


Proof. For A E M we define the operator 


A{X):=A + X\f){f\ (3.23) 

and write a;(A) for the Ith eigenvalue counted from above and ipi{X) for the corre¬ 
sponding eigenvector. Moreover, we remark that a/(l) and </?;(!) correspond to /?; 
and tpi. Assumption (|3.3I) and the definite sign of the perturbation imply that the 
eigenvalues of A(A) are non-degenerate for all A E [0,1]. Thus, standard results, see 
[RS78i Chap. XII], give differentiability of the eigenvalues for all A E (0,1) and we 
apply the Feynman-Hellmann theorem, see e.g. [IZ88] . to deduce for alH E N and 
A E (0,1) 

a\{X) = \{^i{X),ct>)\\ (3.24) 

Hence, we compute using the latter 

(3.25) 



The eigenvalue equation implies 

_ ^ \ 1 

^ - 1 ( 3 . 26 ) 

where we used Fubini’s theorem to interchange the integral with the sum and the 
fact that the vectors {tiW) form an ONB. The resolvent identity (I3.18p implies 


1 -b x{4>, ^ ^ 


A(A) — 2 ; 


(3.27) 
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Therefore, we continue 

eM=- 1 "(■*.( 

J_^2 

iiJ 


1 


( ^)V) d 

Jo 


+ A((/>, -^(p) 

1 

dA \yl + X{(j), -^rzfp) 


1 - 


1 


(</>,^</>) V' Vl+ ^ 

Equation (|3.27p with A = 1 provides the assertion 


1+ ' 


(IT^ = -( 


1 


1 


A-z B-z 


4>)- 


(3.28) 

(3.29) 

□ 


4. Proof of Theorem 12.11 


We decompose the determinant according to the angular momentum decomposition 
|. This implies 

T TT |det((^,qa<«ho)) , 

l^j,k^NL(.E) -!:A I ^VlsCj,fcsCiV[(E) 




(4.1) 

where and correspond to the radial part of the eigenfunctions lying in 

the f-th angular momentum channel and N^{E) to the relative particle number in 
the f-th angular momentum channel. More precisely, 

Ni{E) := #{ /^ G N : 3 j e {1, ■ ■ ■ , Nl} with = Xf } (4.2) 

where (A^(^))^gj,^ denote the eigenvalues of /i£. Since we chose the eigenfunctions of 
Hl and Ha^i to be the same in every angular momentum channel f ^ 1 we obtain 
that only the f = 0 term in the product (14.11) is different from 1. Hence, 


det 


((d>b) 


= det( (<y9:['(0),'!/)^(0))) , (4.3) 

Thus, we reduced our problem to a problem on the half-axis, where the relative 
particle number satisfies 


Lemma 4.1. Given E > 0. Let L and Nl{E) G N such that 


L —>■ oo. Then, 

as L —)• oo. 

Proof. For any E > 0 


NliE) ^ VE 


L 


vr 


=: PoiE), 


p{E) as 

(4,4) 


L^oo |.Br,(0)| L^oo |Hi(0)| 


(4.5) 
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where the first equality follows from e.g. [RS78] . Hence, we obtain for an arbitrary 
e > 0 the inequalities 

Nl{E) ^ #{ A: : ^ i? + e} (4.6) 

for L large enough. Since p is is strictly increasing, we obtain —)■ E. There¬ 

fore, ^^0 (g)(0) E as well because otherwise there would be a gap in the spectrum 

of hP by the definition of the relative particle number N^{E). This implies for an 
arbitrary e > 0 and L large enough 


NliE) 

#{A::A^(0)^i?} 

cp 

ff{k:{^) ^{E-e,E + e)] 

L 

L 


L 



(4.7) 


for some constant c. Since k : A^(0) ^ E}/L —)• po{E), as L —)• oo, this yields 
the claim. □ 


Given (j4.1l) and Lemma 14.1[ Theorem 12.11 will follow from 
Theorem 4.2. Let E > 0. Then, 


det((<y9^(0),V'fc(0)>) 


l^j,k^NL 


= l-C{E)+oil) 


as L ^ oo, Nl G N and N^/L —)• where 

m := ^Sli^/E) 

and 6a is given by Definition \2.2[ 


(4.8) 


(4.9) 


From now we shorten the notation and drop the 0 and L-index of the eigen¬ 
functions and eigenvalues. 

Apart from the product representation discussed in Section [3] the main ingredi¬ 
ent to the proof of Theorem 14.2l is a elementary formula expressing the non-negative 
eigenvalues of the perturbed operator h^ g in terms of the eigenvalues of the opera¬ 


tor /ig plus corrections depending on the scattering phase shift 6a- First, note that 
the eigenvalues of hg can be computed explicitly, see [RS78] . i.e. for n G N 


An = 


/nvrx 2 

\T) 


(4.10) 


Lemma 4.3. Let 6a be given by Definition \2.2[ Then, 

(i) for a ^ 0 and n G N the nth eigenvalues of /ig and g satisfy 

0^ (4.11) 

(ii) for a ^ 0 and n > 1 the nth eigenvalues of /ig and hP^ g satisfy 

(4.12) 
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(iii) and 5 exhibits the following expansion 


= &(v^) - 


+°'i 


(4.13) 


which is valid for all [Xn ^ 0 , and the error term depends on a but is 
independent of n. 


Proof. Let A: > 0. Consider the eigenvalue problem 

— u'l. = kfuk, —47rattfc(0+) + u)j(0+) = 0. 

Introducing Priifer variables 

Uk{x) = Pu{x) sin( 6 »fc(x)) u'^{x) = kpu{x) cos( 6 »fc(x)), 
we see that any non-zero solution of (14.1411 is of the form 

k 


Uk{x) := a sin (^x + arctan 


(4.14) 

(4.15) 

(4.16) 


for some 0 7 ^ a G C. Since any eigenfunction Uk to an eigenvalue k"^ of hf^^^ is a 
solution of (14.1411 in (0, L) and additionally satisfies Uk{L—) = 0, we obtain that 

rife(L) = a sin -|-arctan =0. (4-17) 

On the other hand, all such that (|4.17p is satisfied are eigenvalues of hf^ Since 
the function k e-?• kL + arctan (jfjj) is strictly increasing we obtain for any n G N 
an unique eigenvalue 0 oi h^ ^ such that 


y/Jx^L + arctan 


V dvra ) 


nvr, 


(4.18) 


where pi < p 2 < ''' ■ This proves (i). For the case a < 0 note that ^ admits a 
single negative eigenvalue. Therefore, (|4.18p is only valid starting from the second 
eigenvalue of This implies for all n G N 

_ arct an () _ TT — arctan ( 

- I ^ -(4.19) 

(iii) follows directly from (i), (ii) and Definition (|2.2p from the phase shift. □ 

Corollary 4.4. The eigenvalues of and h^ satisfy 

/ri < Ai < /i2 < A 2 < • • • • (4.20) 

Proof. Note that |(1 q(A:)| < vr for all k > 0. Thus, (14.1011 and (14.1211 imply the 
corollary. □ 

Next we apply the results from Section [3] to the determinant: 

Lemma 4.5. Let G N. Then, 

\pk AjI |Afc pj\ 


det 


N 


l^j,k^N 


n n 

j=l k=N+l 


I ^k I I pk pj I 


(4.21) 
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Proof. First note that ^ is bounded from below by ( 12 . 71 ) . This and ^ 0 
imply —E^ p{hL) r\ p{h^ l) foi" some E > 0. Moreover, ( 12 . 71 ) provides 

1 1 


E.a I 


hr + E 


hr T + E 




(4.22) 


for some E L^((0, L)) and Corollary 14.41 gives 


1 


> 


1 


> 


1 


> 


1 


Pi + E Ai + i? p2 + E X 2 + E 


> 


(4.23) 


W+E 


and 


the eigenvalues satisfy assumption O. Furthermore, the operators 
■ 0 ^ g are non-negative with trivial kernel and compact. Therefore, we are in po- 
sition to apply Theorem 13.11 and obtain 




N 


l^j^k^N 


n n 


1 1 

1 1 

1 1 

1 1 

1 /Afc+^ 

\j+E \ 

1 

Hj -\-E \ 

1 1 

1 1 

1 1 

1 1 

1 

\j+E \ 

1 (J-k+E 



N 


n n 

j=i ^=^- 1-1 


\Pk Pj\ 

I Pk Pj I I Pk pj I 


(4.24) 

□ 


Proof of Theorem 14.21 We start with the product representation given in 
Lemma 14.51 Note that for a < 0 there is an ambiguity since there exists pre¬ 
cisely one negative eigenvalue pi. Therefore, we treat the j = 1 term in the product 
separately. We define 


•4? = n 


\hk — All |Afc — /III 


n h 


ihk — Aa:)(Ai — pi) 


k=N+l 1'^*' hl\ k=N+l 


and estimate using Corollary 14.41 
{pk Afc)(Ai /ii) 


E 

fc=7V-|-l 


(Afc — Xi){pk — pi) 


^ iAi — pi\ 

k=N+l [[ — ) - [l) 


(Afc — Xi){pk — pi) 

((^)2_ (11^)2) 


( 4 , 25 ) 




{2k - 1) 




< c 




(4.26) 


Since is uniformly bounded from below with respect to L, see Lemma 12.9 


InM^ 



\pk All |Afc ^i| A _ 

\\k-\i\\Pk-pi\)~ 


(4.27) 


as A^, L —>■ 00 and ^ —)• p{E) > 0. Therefore, we are left with a product consisting of 
the non-negative eigenvalues and apply Lemma 14.51 use Lemma 14.31 (i) and y/XP = 
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n G N, to obtain 


In 


l^j,k^N 


= lnA^ 


+ V V in(^ l \\ iikT^yf - 0'^ - \ \ ,^28) 

j=2k=N+i ^ |(^7r)^-(j7r)^|| (A:7r-(5(v/7I^))^-(j7r-J(0I-))^| J 

In the following the 0(1) and o(l) terms refer to the asymptotics L,N ^ oo, 
N/L —>• po{E) > 0. Equation (I4.27p above, Lemma [A. II below and the abbreviation 
Qk ■= for /c G N yield 

rt™ = _V V (2jgj+g|) (2fcgfc + gD 

kh'+i (+ 3kf - U + Qj? ) - j2) 

Using Lemma [A.2I and the abbreviation 6k ■= —^6{^/Xk) for /c G N, we have 

^ °° {2j6j + 6]) {2k6k + 6l) 


(4.29) 


irai) = - V V -— ' V ' — o —- + o(i). 

Lemma IA.3I implies 

N 2N , r r 

H 33 = -E E 7 ^^ + oW- 

j=2 k=N+l - J > 


Lemma IA.4I yields 

iV 

(I4.3ip = —\ [ dx [ 

Jo Jt 


‘ixy6a{x'x)6a{yT^) , 
dy———1717-^ 0(1)- 


We define for 0 ^ x < y 


g{x,y) := 


iV+l (y 2 — ^ 2)2 

4xy<5a(7rx)(I„(7ri/) 


{y + x)2 

The explicit representation of 6a implies for all e > 0 

sup sup ||(V 5 r)(x,i /)||2 := c(e) < oo. 

6>e (3:,y)G(0,6)x (fo,oo) 

Therefore, using the mean value theorem and the Cauchy-Schwarz inequality, we 
compute for a 0 < e < y/E and iV, L big enough 


(4.30) 

(4.31) 

(4.32) 

(4.33) 

(4.34) 




2 N 

L 


dy 


L 

N ^ 2 N 


‘ixy6a{x'K)6a{yT^) ^2 


-Si (N/L) 


1 


(y + x)2 “ {y- x)2 

1 


^c(e) [ dx [ dy||(A^/L-x,y-iV/L)| 

Jo J^ + i- 


^2c(e) 


iV 

L 


dx 


'0 


L 

r 

L 


'{y - xf 


dy 


{y - x) 


= 0 ( 1 ), 


where we used the inequality 

|x-iV/L| + |y-lV/L| 


(y - xy 


^ 2 


{y-xY 


(4.35) 


(4.36) 
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which is valid for all x < N/L < y. Moreover, since ^ —)• > 0, we compute 

- = lnL + 0(l). (4.37) 

) 

Hence, combining equation (14.351) and (|4.37p . we end up with 



g32D =-InL^r^^CvriV/L)+ 0(1) (4.38) 

TT^ 

=-lnL\6l{VE)+oilnL), (4.39) 

where the last line follows from tt^ —>■ ^/E. This gives the assertion. □ 


Appendix A. Proof of the auxiliary lemmata 


In this section we prove the missing lemmata used in the proof of Theorem 14.21 
We do not claim to give optimal or very elegant estimates. Throughout this section 
we drop the index a in the scattering phase shift and restrict ourselves to the case 
a < 0. This implies the following estimate on the phase shift 

(5(x) - (5(?/) ^ 0, (A.l) 


for X < y, which we use in the sequel. The case a ^ 0 is even simpler since in that 
case the Definition (12.2p of the phase shift implies the uniform bound 

ll<5||oo ^ (A.2) 

which simplifies some of the following estimates. Moreover, we use the elementary 
asymptotics 


N CO 

E E 

j=l k=N-\-l 
N CO 


(k-jy 


E E 

j=l k=N+l 


{k - j)f^ 


= O(lniV), 


= 0 ( 1 ) 


as A —>■ oo, where /3 > 2. 

Lemma A.l. Set := —k5{yfJjLk) for /c G N. Then, 

^ k^+1 \{ik + 9kf - U + 9jf ) (k^ - P) ) 


(A.3) 

(A.4) 


(A.5) 

(A.6) 


as N,L ^ oo, 


TT 
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Proof. We prove the assertion in two steps. First we consider the j = N and 
k = N + 1 summand. Note that Lemma o above and E > 0 imply 


lim gN = ^ lim g^+i = 

N^L —^oo NjL —^cxD 

N/L^^/E/tt AT/L^Ve/tt 


> - 1 . 


TT 


(A.7) 


Thus, for j = N and fc = + 1 

lim In 


/ ((iV + 1 + gN+iV - N^) {{N + 1)^ - (iV + ) 

“ V ( (iV + 1 + 5jv+l)' - (iV + gNf ) {{N + 1)2 - iV2) 

N/L^VE/tt 


= lim In 
N,L^oo 

NIL^y/EI-ir 


(l + ^Af+i) (l — On) (2iV + 1 + gN+i) (2A^ + 1 + gN) 
(l + BN+i — gN) (2A^ + 1 + ^Af+i + gN) (2A^ + 1) 


= In 1 - 


sXVe)\ 


7T^ 




Moreover, along the same line using (IA.7P 

{2NgN + g%) (2(iV + l)5(Ar+i + ^at+i) 


lim — - 


(A.8) 

SHVE) 


nX^^ ((iV + 1 + gN+if -{N + gN? ) {{N + 1)2 - N^) 

NIL^y/EI-K ^ \ / 


Therefore, the j = N and k = N + 1 term is of order 1. 
For j^A^<A^+l<A:we want to apply the bound 

.2 1 


TT^ 


ln(l + x) — X ^ 


2 1 - X 


for X G M with |x| < 1, to x = Xjk where 

(2jffj + g]) {2kgk + gl) 


^jk • — 


{{k + gk?-{j + gj)‘'){k^-P) 
We estimate using \gn\ ^ 1 for all n G N and gi^ — gj ^ 0 

I k I ^ 


(A.9) 


(A.IO) 


(ATI) 


(2i + gj){‘2k + gj.) 

1 

(j +gj + k + gk){k + j) 

{k- j + gk -gj){k- j) 


< (A,12) 

Since j ^N<N+l<k, this implies in particular \xjk\ ^ and we continue 
using (lA.lOp and (IA.12I1 

N OO N OO 

E E I ln(l + Xjk) - Xjk\ <E E Xjk 

j=l k=N+2 j=l k=N+2 

N OO T A 

<E E ijX) - 0 ( 1 ). 


j=2 k=N+l 

as N ^ OO, where we used (IA.4h in the last line. 


(A.13) 

□ 
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Lemma A. 2. Define 6k ■= —^S{y/Xk) for fc G N. Then, 
^ °° {2j6j + 5 |) {2k6k + 6l) 


E E 

j=2 k=N+l 


as N,L ^ oo, ^ —>■ 


((A: + 6k)‘^ — {j + 6j)‘^ ) 
{2jgj + g]) {2kgk + gl) 


({k + gk) - {j + gj) ) 


{k - jf 


= 0 ( 1 ) 


(A.14) 


Proof. First, using the expansion of Lemma [4.31 we obtain for all n G N, n > 1, 

(A.15) 


bn-^nl ^ ^ ll-^llooP'l 


ttL 


c 

V 


where the constant c > 0 depends only on a. We prove the assertion in two steps. 
In the first step we consider the numerator only in the second step we consider the 
denominator. Using (jA.lhp we estimate 

^ °° ' + ^]) {‘^kSk + 6“^) — {2jgj + gj) {“^kgk + g"^ 


E E 

j=2 k=N+l 
^ N oo 

^-E E 

j=2 k=N+l 
N oo 


{{k + gk?-{j+g,f){k^-P) 
{j + l){k + l) 




L 

C 


{ik + gkf-{j + gjf){k^-P) 


tE E 


(j + l)(fc + l) 


j=2 k=N+i (^ + J - 2) {k + j) {k - jf 

as N,L ^ oo, ^ ^ where we used \gj + gk\ ^ 2, gk — gj >0 for j < k and 
(]A.3I) . In order to estimate the denominator we use (lA.lSP to obtain some constant 
c > 0 independent of j, k such that 

k + j 


= O 




(A.16) 


((A: + gkf - (j + gjf'j - ((A: + 6kf - {j + 6jf^ 


< c- 


L 


(A.17) 


Thus, 


N OO 

{2j6j+6]){2k6k + 6l) 

j=2 k=N-\-l 

1 


{ik + gkf-{j + gjf){k^-jf 


. N oo 

s-E E 


((A; + 6kf - (j + 6jf) [kfi - jf 


jk{k + j) 


L 


j^ 2 kfN+i {k‘^-jf {{k + gkf - {j + gjf) {{k + 6kf - {j + 6j)f 

\T 

jk 


. N OO 

^^E E 


L 


'l^2kfN+i {k-jf{k + j-2Y{k + j) 


3 


= 0 ( 1 ) 


(A.18) 


as N,L ^ oo N/L where we used I^a, + gj\ ^ 2, \6k + ^ 2, gk — > 0 

and 6k — 6j > 0 for j < k. □ 
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Lemma A.3. The estimate 

N CO 


E E 

j=2 k=N+l 


(■ijSi+S]) (2kS,+Sl) 

{(fe + 4)^ - (i + 4)^) T - f^) 


N 2N 

E E 

j=2 k=N+l 


Ajk5j5k 




^ = 0 ( 1 ) 

(A.19) 


holds as N,L ^ oo, p 

■ ’ I, TT 


Proof. First, we bound the tail, i.e. using 6^ — 6j > 0 for k > j and |(5ri| ^ 1 for 
all n E N we estimate 

N OO (oAs: I (nUK \ ^ oo 


E E 


j=2 k=2N+l 
oo 

^ ^ (k-NY 

k=2N+l ^ ’ 


{2j5j + 6]) { 2 k 6 k + si) 

({k + 6kf - {j + djf ) (/c^ - P) 


€ 


E E 

j=2k=2N+l 


1 


{k-jf 


N 


= 0 ( 1 ), 


^jk5j5k 


N 2N 

as A —>■ OO. We insert zb > > „ „x 

U fcir+1 ((^+- j ) 

in the next step 6k — > 0 yields 

^ 2iv (2j<5,- + 6]) {2k6k + 51)- Ajk5j 4 


(A.20) 
in (lA.lOp . Thus, 


E E 

j=2 k=N+l 
N 2N 

<E E 

j=2 k=N+l 
N 2N 

<3E E 

j=2 k=N+l 


{{k + 6kf-{j + 6jf){k^-p) 
2{k + j) + l 


[k-j) {k + j){k + j - 2)) 

' =0 


{k-j)\k + j-2) 

as A —>■ oo, where we used (IA.3P in the last line. In the third step, again |(5n| ^ 1 
for n E N yields 

N 2N 

E E 


In A 

~ir 


(A.21) 


Ajk 


j=2 k=N+l 
N 2N 

<E E 


(fc2 - p) 


1 


1 


{{k + — {j + 6jY) (A;2 P) 

9jk (k + j) 


J 


i=2 k=N+i - P) {k + j- 2) {k - j) 


N 2N 

<sE E 


1 


j=2 k=N+l {k j) 
as A —7> OO, where we used ()A.4h . 

Lemma A.4. The asymptotics 
N 2N 


= 0 ( 1 ), 


E E 


Ajk6j6k 


2 - i ^ dx 


' 2N 
L 


dy 


4xy(5(x7r)(5(y7r) 


holds as N,L ^ oo, p —?• 

’ ’ L n 


(A.22) 

□ 

= 0(1) (A.23) 
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Proof. We recall that 6k ■= —-(5(\/^) and we rewrite 


AjkSjdk 


N 2N 

y y 

j=2 k=N+i y ~ y 

Thus, we estimate 


N 2N 

E E 




12^2/.^ U _ I 

j=2k=N+l \\l) \l) ) 


(A.24) 


1 


N 2N j kxfj 


E E 


L L 


<5(¥)hx) 


N 2N 


//"" L 


2N+1 


L xy5ixTT)5iyTT) 
dy 


^ {y^-x^y 


k + l 


j=2k=N+l L 


where 


y) xy6{x7r)6{y7r) 


(A.25) 


(?/2 - xy 

Using the mean-value theorem and the Cauchy-Schwarz inequality we obtain 
N 2N 

^y y sup |(V/)(x,y )|2 

j=2 k=N+l U-S/)6(^,i)x(i^) 


(A.26) 


dx 


fc+i 

L 


f .?-l 

L 


dy 


L-^^L-y 


N 2N 




sup 


^ j=2k=N+l{x,y)&{i^,j-)x{^,!^) 

where 1-12 denotes the Euclidean norm. We compute 

1 


|(V/)(x,i/)| 


2 ’ 


(V/)(x,2/) = 


(y2 - xy 

(y^ — x‘^){y6{x7r)6{yTr) + xyS' {xTT)6{yTT)Tr) + 4:x‘^y6{x7r)6{y7r) 
(y^ — x‘^){x6{xTr)6{yTr) -|- xy6{xTT)6'{y'K)TT) — 4:xy‘^6{x7r)6{yTr) 

1 

■9ix,y)- 


(A.27) 


(A.28) 


(y2 _ j.2^^ 

We estimate for (x,y) E i) x (xi j ^ ^ 

1 


y^ — 




(A:-bj-l) {k-j) 


< 




L6 


_ -y 


(A.29) 


(A.30) 


and, using 6,6' € L°°((0, oo)). 


sup 


15(2:,y )|2 ^ 


sup 

N+ 

L 


(x,j/)e(o,™)x(o,^) 


|5'(3:>y)|2 = 0(1) (A.31) 


as N,L ^ 00 , ^ Thus, ()A.30I1 and ()A.3ip imply 




N 2N 


j=2 k=N+l 


(y - x)" 


= 0 ( 1 ) 


(A.32) 
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as N,L ^ oo, -jj 


Ve 

TT 


□ 
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